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ABSTRACT 

Context.  Eruption  of  a  coronal  mass  ejection  (CME)  drags  and  “opens”  the  coronal 
magnetic  field,  presumably  leading  to  the  formation  of  a  large-scale  current  sheet  and 
the  field  relaxation  by  magnetic  reconnection. 

Aims.  We  analyze  physical  characteristics  of  ray-like  coronal  features  formed  in  the 
aftermath  of  CMEs,  to  check  if  the  interpretation  of  this  phenomenon  in  terms  of 
reconnecting  current  sheet  is  consistent  with  the  observations. 

Methods.  The  study  is  focused  on  measurements  of  the  ray  width,  density  excess,  and 
coronal  velocity  field  as  a  function  of  the  radial  distance. 

Results.  The  morphology  of  rays  indicates  that  they  occur  as  a  consequence  of  Petschek- 
like  reconnection  in  the  large  scale  current  sheet  formed  in  the  wake  of  CME.  The 
hypothesis  is  supported  by  the  flow  pattern,  often  showing  outflows  along  the  ray, 
and  sometimes  also  inflows  into  the  ray.  The  inferred  inflow  velocities  range  from  3  to 
30  kms-1,  consistent  with  the  narrow  opening-angle  of  rays,  adding  up  to  a  few  degrees. 

The  density  of  rays  is  an  order  of  magnitude  larger  than  in  the  ambient  corona.  The 
density-excess  measurements  are  compared  with  the  results  of  the  analytical  model  in 
which  the  Petschek-like  reconnection  geometry  is  applied  to  the  vertical  current  sheet, 
taking  into  account  the  decrease  of  the  external  coronal  density  and  magnetic  field 
with  height. 

Conclusions.  The  model  results  are  consistent  with  the  observations,  revealing  that  the 
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main  cause  of  the  density  excess  in  rays  is  a  transport  of  the  dense  plasma  from  lower 
to  larger  heights  by  the  reconnection  outflow. 

Key  words.  Sun:  coronal  mass  ejections  (CMEs)  -  Sun:  corona  -  (Sun:)  solar  wind  - 
magnetohydrodynamics  (MHD) 


1.  Introduction 

According  to  the  current  comprehension  of  solar  coronal  mass  ejections  (CMEs),  the  erup¬ 
tion  of  an  unstable  magnetic  structure  is  tightly  associated  with  the  formation  of  the  large 
scale  current  sheet  (hereinafter  CS)  in  the  wake  of  the  eruption.  This  concept,  connect¬ 
ing  the  large  scale  eruption  and  the  localized  energy  release  in  a  form  of  flare,  was  put 
forward  and  initially  developed  by  Carmichael  1964;  Sturrockil966t  Hiravamall974t  Kopp 
&  Pneuman  1 19761  (the  so-called  CSHKP  model).  When  the  current  sheet  becomes  long 


enough,  the  tearing  instability  sets  in,  leading  to  fast  reconnection  of  magnetic  field  (for 
the  analytical,  numerical,  laboratory,  and  observational  results  see,  e.g., 


Ugai 

1987; 

Gekelman  &  Pfister 

1988 

Furth  et  al 


1963 


Vrsnak  et  al 


20031.  respectively).  Reconnection  in 


the  post-CME  current  sheet  results  in  an  abrupt  energy  release  causing  a  flare,  and  on  the 
other  hand,  enhances  and  prolongs  the  acceleration  of  the  erupting  magnetic  field  struc¬ 
ture  (e.g., 


Lin 


2004; 


Vrsnak 


2008,  and  references  therein).  Numerical  MHD  simulations 


covering  various  scales,  from  low  corona  up  to  1  AU,  also  demonstrated  post-CME  cur- 


Torok  et  al. 

2004; 

Rilev  et  al. 

2007 

Rilev  et  al. 

2002: 

Roussev  et  al. 

2003 

20071,  for  a  review  see  Forbes  et  al.  2006). 


The  most  prominent  consequence  of  CME-associated  reconnection  is  the  appearance  of 
the  so  called  two-ribbon  flare,  closely  synchronized  with  the  CME  acceleration  stage  (e.g., 


Maricic  et  al. 


20071  and  references  therein).  The  ribbon  expansion  away  from  the  magnetic 


inversion  line,  associated  with  growing  system  of  hot  loops  connecting  the  ribbons,  has  led 
to  the  formulation  of  the  CSHKP  model.  The  discovery  of  cusped  structure  above  the  flare 

looptop  hard  X-ray 


loops  (e.g., 
sources  (e.g., 


Tsuneta  et  al 


1992; 


Masuda  et  al 


Forbes  fc  Acton 


2004; 


1994; 


2006),  the  loop  shrinkage  (e.g., 


Vrsnak  et  al 


1996 


Aschwanden  et  al 


Svestka  et  al 


1987 


Tsuneta 


1996; 


1996) 


Sui  et  al 


Forbes  fc  Acton 


2004 


1996; 


Verpnig  et  al 


Sheelev  et al 


20061).  and  the  recognition  of  growing  posteruption  loop  systems  in  the 


absence  of  flares,  further  supported  the  CSHKP  scenario.  Additional  evidence  for  recon- 


(e.g.,  Webb  &  Cliver  1995; 

simnett  et  al.  1997: 

Wang  et  al. 

1999) 

downflows  above  the 

post-eruption  area 

horizontal  converg 

des  (e.g.,  I 

TcKenzie  &  Hudson 

1991 

:  Innes  et  al.  2 

!003:  Asai  et  al. 

2004), 

ing  flows  above  the  loops  (Yokovama  et  al 

2001 

Lin  et  a' 

.  2005 

),  and 

flare-associated  radio  emission  below  the  eruptive  prominence 

(Vrsnak  et  al. 

>003). 

Recently,  post-CME  features  seen  in  UV  spectra 

Uiaravella  et  a 

1.  2002; 

Ko  et  al 

2003; 

Innes  et  al. 

2003b 

Lin  et  al 

2005;  Bemoorad  et  al. 

2006;  C 

laravella  &  Ravmoncl 

2008), 

X-rav  images  (Sui 

et  al.  200z 

1.  2005),  and  white  light  coronagraph  images 

Ko  et  al 

2003; 

Webb  et  al. 

2003; 

Lin  et  al. 

2005)  have  been  attributed  to  the  current  sheets  expected  to 
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“connect”  the  flare  loops  to  the  CME  core.  The  UV  spectral  signature  is  generally  emis¬ 
sion  in  the  high  temperature  lines  of  Fe  XVIII  or  FeXIX  formed  at  5-  10  MK.  EUV  images 
from  the  Extreme-ultraviolet  Imaging  Telescope  (EIT)  on  board  the  Solar  and  Heliospheric 
Observatory  (SoHO)  and  the  Transition  Region  and  Coronal  Explorer  (TRACE)  show  emis- 


sion  in  FeXXIV  at  even  higher  temperatures  ( McKenzie  &  Hudso 


1 


The  UV  features  decay  slowly,  on  time  scales  from  several  hours  (|Ciaravella  fc  Raymond 


200811  up  to  a  few  days  (Bemporad  et  al 


2006 


Ko  et  al. 


993; 


Innes  et  al.  20031 


2003).  The  X-ray  data  from  the 


Reuven  Ramaty  High  Energy  Solar  Spectroscopic  Imager  (RHESSI)  show  regions  of  high 
temperature  emission  above  loop  tops,  and  both  the  morphology  at  different  energies  and 
the  evolution  of  the  structures  supports  the  CS  indentification.  The  white  light  features 
are  bright,  narrow  rays  that  appear  to  map  from  the  cusped  flaring  loop  system  to  the 


CME  core  (Lin  et  al 


2(105).  Some  show  blobs  of  plasma  moving  outwards  at  several  hun¬ 


dred  kms-1.  Like  the  UV  features,  the  white  light  structures  can  last  for  a  day  or  two.  In 
general,  white  light  observations  reveal  the  electron  column  density,  UV  data  show  tem¬ 
peratures,  emission  measures  and  Doppler  shifts,  and  X-ray  observations  reveal  very  high 
temperature  plasma. 

The  aim  of  this  paper  is  to  demonstrate  that  post-CME  rays  appear  as  a  consequence 
of  the  reconnection  in  the  current  sheet  formed  in  the  wake  of  CME,  and  to  quantify  their 
basic  physical  characteristics.  We  outline  the  working  hypothesis  in  Sect.  2,  whereas  details 
of  the  model  used  to  compare  the  theoretical  expectations  with  observations  are  provided 
in  the  Appendix.  Observations  and  measurements  are  presented  in  Sect.  3,  focusing  on 
the  morphology  and  plasma  densities  in  post-CME  rays.  The  temperature  structure  and 
the  ionization  state  will  be  presented  in  a  separate  paper.  Our  density  measurements  are 
compared  with  the  model  estimates  and  some  previous  empirical  results  in  Sect.  4.  Results 
are  discussed  and  conclusions  drawn  in  Sect.  5. 


2.  Working  hypothesis 

Since  post-CME  rays  extend  outward  from  cusped  structures  associated  with  growing  post¬ 
eruption  loop-systems,  our  working  hypothesis  is  that  they  appear  as  a  consequence  of 
the  reconnection  in  the  CS  formed  in  the  wake  of  CME.  Furthermore,  we  assume  the 
reconnection  process  can  be  approximately  described  in  terms  of  the  steady-state  Petschek 
regime  ( PetscheiJl9641.  In  this  regime,  the  reconnection  takes  place  within  a  small  diffusion 
region  (DR),  out  of  which  two  pairs  of  standing  slow- mode  shocks  (SMSs)  extend  along 
of  the  axis  of  symmetry.  Thus,  the  CS  is  bifurcated,  with  SMSs  separating  the  inflow  and 


outflow  region  ([Petschek 


1964: 


Soward  &;  Priest 


1982).  At  SMSs  the  inflowing  plasma  is 


compressed,  heated,  and  accelerated,  forming  the  upward  and  downward  reconnection  jet 
(Fig.  CD).  Note  that  electric  currents  are  concentrated  only  in  DR  and  SMSs. 

When  reconnection  takes  place  in  the  homogeneous  environment,  the  characteristics  of 
outflowing  plasma  are  determined  primarily  by  the  external  plasma-to-magnetic  pressure 


Appendix  in 


Aurass  et  al. 

2002: 

Skender  et  al. 

CO 

O 

O 

CM 

200311.  However,  in  the  case  of  a  vertical 


current  sheet,  the  ambient  coronal  density  and  magnetic  field  decrease  with  height,  so 
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A1 

I 


% 


Fig.  1.  a)  Interpretation  of  the  ray  in  terms  of  the  Petschek  (|1964f)  reconnection  model: 
DR  -  diffusion  region;  SMS  -  slow  mode  shocks  (dashed  lines  intersecting  at  DR);  PEL 
post-eruption  loops.  Magnetic  field-lines  are  drawn  by  thin  lines;  gray  area  between  the 
SMSs  outlines  the  post-CME  ray.  The  coordinate  system  is  indicated  (the  line  of  sight  is 
in  y-direction).  b)  An  element  of  bifurcated  current  sheet;  arrows  mark  plasma  flows  (uin 
inflow  velocity;  vsw  -  solar  wind  speed;  uout  -  outflow  speed),  where  the  arrow  thickness 
depicts  plasma  densities. 


Fig.  2.  Dependence  of  the  width  and  morphology  of  the  ray  at  a  given  height  on  the  ori¬ 
entation/geometry  of  the  current  sheet:  a)  straight  CS;  b)  wavy  CS.  The  CS  half-thickness 
is  denoted  by  d,  horizontal  length  (i.e.,  the  length  in  zy-plane,  see  also  Fig.  1)  by  L,  the 
ray  width  by  D,  and  the  inclination  to  the  line  of  sight  (LoS)  by  6.  The  plasma-column 
length  A  is  drawn  by  thick  double-arrow.  At  the  bottom  the  intensity  profile  is  sketched. 


the  characteristics  of  the  reconnection  jets  depend  also  on  the  height.  The  main  effect  is 
transport  of  the  dense  plasma  from  lower  heights  upward,  making  the  outflow  jet  much 
denser  than  the  ambient  corona  (see  Appendix).  In  this  respect,  the  post-CME  ray  should 
be,  to  a  certain  degree,  similar  to  coronal  streamers. 
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Bearing  in  mind  the  geometry  of  the  Petschek-type  reconnection,  the  post-CME  CS 
should  show  distinct  morphological  characteristics.  It  should  be  thinnest  at  the  height 
of  the  diffusion  region,  getting  wider  with  increasing  height.  In  the  homogeneous  plasma 
the  angle  between  the  SMSs  is  determined  by  the  inflow  Mach  number  llSoward  fe  Priest 


982; 


Vrsnak  fe  Skendei 


(Vrsnak  fe  Skender 


2005).  the  half-angle  typically  adding  up  to  several  degrees 


20051).  Thus,  for  the  CS  length  of  several  solar  radii,  the  Petschek-CS 


thickness  (the  distance  between  the  SMSs)  is  on  the  order  of  100  Mm.  Note  that  in  the  case 
of  the  vertical  CS  in  the  solar  corona,  where  the  magnetic  field  diverges  radially,  the  angle 
between  SMSs  should  increase  with  the  height  even  at  uniform  inflow  Mach  number  (see 
Appendix),  so  the  CS  thickness  should  be  somewhat  larger  than  in  the  plane-symmetric 
case. 

Further  important  characteristics  of  the  post-CME  CS  that  determine  how  it  should 
look  in  coronagraphic  images,  are  the  horizontal  length  and  the  orientation  with  respect  to 
the  line  of  sight  (Fig.  [2]).  In  this  respect,  it  is  important  to  note  that  the  angle  9  between 
the  plane  of  CS  and  the  line  of  sight  (LoS)  has  to  be  rather  small,  otherwise  the  CS  would 
not  appear  as  a  distinct  feature  in  coronagraphic  images.  This  causes  a  “selection  effect” , 
which  can  explain  the  fact  that  post-CME  rays  can  be  recognized  only  in  a  relatively  small 


fraction  of  events  (Webb  et  al 


2003). 


In  Fig.  2a  we  depict  how  the  effective  thickness  D  of  the  ray  depends  on  the  CS  thickness 
2d,  horizontal  length  L,  and  the  angle  9.  On  average,  a  ray  should  be  thinner  and  brighter 
for  smaller  9  (larger  LoS  column-length  A).  In  Fig.  [5}j  we  illustrate  how  the  variation  of 
9  along  the  LoS  can  lead  to  a  complex  appearance  of  the  ray  (see  also  Saez  et  al.  2007). 


i.e.,  a  wavy  CS  would  be  seen  by  observer  as  a  “multi-ray”  structure.  Note  that  the  angle 
9,  or  its  variation  along  the  LoS,  can  vary  in  time  and  height,  so  the  morphology  of  the 
structure  can  change  in  time  and  can  be  different  at  different  heights. 


3.  Observations 


3.1.  SoHO/LASCO  mass  images 


For  the  analysis  of  the  morphology  and  the  column-density  of  post-CME  rays  we  employ 
data  provided  by  the  Large  Angle  Spectroscopic  Coronagraph  (LASCO; 


Brueckner  et  al 


19951)  on  board  SoHO.  We  focus  on  the  LASCO-C2  images,  covering  the  radial  distance 
range  from  2  to  6  solar  radii  (r®).  In  particular,  we  use  the  so-called  mass  images  that  are 
derived  from  LASCO  base-difference  images  calibrated  in  units  of  solar  brightness.  These 
images  show  changes  of  the  column-mass  along  the  line  of  sight  under  the  assumption  that 


the  structure  lies  in  the  plane  of  sky  (for  details  we  refer  to  iBillingsl Il96f ; 


Poland  et  al 


1981 


Vourlidas  et  al 


200C 


20021) ;  bright  pixels  (positive  values)  represent  areas  where  the 


column-mass  is  increased,  whereas  dark  pixels  (negative  values)  show  depleted  regions. 

The  initial  sample  of  post-CME  rays  considered  for  the  analysis  consisted  of  eleven 
events.  We  emphasize  that  these  rays  were  spotted  while  analyzing  the  associated  CMEs 
for  other  purposes,  i.e.,  the  rays  were  chosen  rather  randomly.  Thus,  no  attempt  was  made 
to  analyze  the  occurrence  rate  of  post-CME  rays  (for  the  criteria  and  statistical  background 
we  refer  to 


Webb  et  al 


2003).  Nevertheless,  it  should  be  noted  that  ray-like  features  in  the 
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wake  of  CMEs  are  not  uncommon,  particularly  bearing  in  mind  specific  circumstances 
under  which  a  post-CME  current  sheet  could  be  observed.  That  is,  the  level  of  post-CME 
activity  should  not  be  too  high,  current  sheet  has  to  be  oriented  at  small  angle  with  respect 
to  the  line  of  sight,  it  should  be  relatively  stable,  etc. 

Inspecting  the  selected  sample  we  noticed  intermittent  activity  in  most  of  the  rays, 
manifested  as  changes  of  their  shape,  contrast,  inclination,  etc.  These  changes  were  either 
caused  by  coronal  disturbances  coming  from  remote  eruptions,  or  were  revealing  internal 
activity,  such  as  blob  formation/ejection  and/or  outward  propagating  wave-like  perturba¬ 
tions. 

For  a  detailed  analysis  we  chose  three  post-CME  rays,  observed  on  8-9  January  2002, 
18  November  2003,  and  26  June  2005.  The  events  of  18  November  2003  and  26  June  2005 
were  selected  as  examples  of  relatively  stable  rays.  In  both  events  we  analyzed  only  one  ray 
image  (12:50  UT  and  06:30  UT,  respectively).  On  the  other  hand,  the  event  of  8-9  January 
2002  was  showing  a  significant  activity,  similar  to  that  in  the  remaining  eight  events.  This 
event  was  selected  to  quantify  the  level  of  variability  in  post-CME  rays,  so  we  performed 
measurements  at  three  different  times  (00:06,  06:06,  and  12:06  UT).  Note  that  rays  in  the 
selected  exposures  did  not  show  significant  blob-like  features  or  similar  inhomogeneities, 
although  they  might  have  been  present  in  earlier  and/or  later  times. 


3.2.  The  ray  morphology 

In  Fig.  [3]  we  illustrate  the  variability  of  the  ray  morphology  by  showing  two  images  of  the 
ray  that  became  recognizable  a  few  hours  after  the  east-limb  CME  eruption  on  8  January 
2002  (first  appearance  in  LASCO-C2  at  17:54  UT).  In  Fig.  [3jr  we  first  show  a  simple  ray 
pattern  that  was  recorded  half  a  day  after  the  eruption.  Four  hours  earlier,  the  ray  had  a 
complex  form  consisting  of  three  radial  substructures  (Fig.  [3jo);  this  “multi-ray”  pattern 
could  be  a  result  of  a  wavy  current  sheet,  as  illustrated  in  Fig.  [2)j.  Morphological  changes  of 
this  kind,  developing  on  the  time  scale  of  hours  are  a  common  characteristic  of  post-CME 
rays,  since  they  were  present  in  nine  out  of  eleven  events  from  our  initial  sample. 

The  estimates  of  the  ray  width  and  the  column-mass  excess  are  based  on  the  LASCO- 
C2  mass-images  (Fig.  [3] left),  which  provide  for  each  pixel  the  difference  Am  between  the 
column-mass  in  the  actual  image  and  the  reference  image.  In  the  right-hand  panels  of  Fig. [3] 
we  show  two  mass-difference  profiles  A m(y),  where  Am  is  expressed  in  g/pixel,  while  the 
abscissa,  corresponding  to  the  y  coordinate  of  LASCO  images,  is  represented  in  pixels.  The 
profiles  shown  were  measured  in  the  LASCO-C2  difference-images  along  the  bin  marked  by 
white  lines  in  the  left  panel  of  Fig.  [3j  i.e. ,  roughly  perpendicular  to  the  ray  direction.  The 
bin  is  10  pixels  wide,  and  the  profile  shows  A? n  averaged  over  the  bin  width.  Pixels  with 
Am  =  0  represent  areas  where  the  column-mass  remained  the  same  as  in  the  reference 
image.  Note  that  the  mass  astride  the  rays  is  depleted  (Am  <  0),  probably  due  to  the 
CME-associated  coronal  expansion,  usually  seen  as  coronal  dimming  when  the  eruption  is 
launched  from  regions  on  the  solar  disc.  The  depletion  might  be  partly  caused  also  by  the 
presumed  inflow  into  the  current  sheet. 
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Fig.  3.  a)  Simple  post-CME  ray  (9  January  2002  at  06:06  UT);  b)  Complex  ray  consisting 
of  three  radial  substructures  (9  January  2002  at  02:30  UT).  The  ray  is  marked  by  white 
arrow.  The  scale  of  the  LASCO-C2  images  is  given  in  pixels  (1  pixel  corresponds  to  11.9 
arcsec).  The  column- mass  profiles  measured  along  the  bin  marked  by  white  lines  are  shown 
in  the  right-hand  panels  (abscissa  scale  in  pixels).  Estimates  of  the  ray  width,  D,  and 
column-mass  excess,  8m,  are  indicated  in  a).  Dashed  arrows  indicate  the  lower  and  upper 
limit  of  D. 


Given  the  working  hypothesis,  where  the  ray  is  considered  to  be  the  reconnection  out¬ 
flow  jet  of  a  large-scale  Petschek-like  bifurcated  current  sheet,  we  are  interested  in  the 
column-mass  excess  of  the  ray,  8m,  defined  as  the  difference  between  the  peak  in  the  ray 
column-mass  Aror  and  the  column-mass  in  the  adjacent  corona  A mc.  In  other  words,  8m 
determines  the  difference  of  column-mass  associated  with  the  current  sheet  and  the  inflow 
region  (8m  =  A mr  —  Amc;  see  the  right-hand  panel  of  Fig.  [3js.) .  The  ray  width  D  is  mea¬ 
sured  as  the  full  width  of  the  ray  structure  in  the  mass-difference  profile  (Fig.[3ji  right).  The 
main  source  of  error  in  estimating  8m  and  D  is  in  determining  the  value  A mc,  since  this 
value  is  frequently  different  on  opposite  sides  of  the  ray  (Fig.[3ji  right).  Thus,  we  estimated 
the  upper  and  lower  limit  (dashed  arrows  in  Fig.  [3ja-  right),  and  these  measurements  gave 
us  the  error  bars  plotted  in  Fig.  [4] 

The  measurements  illustrated  in  Fig.  [3]  were  performed  in  five  images  of  post-CME 
rays  over  the  LASCO-C2  range  of  heights  to  estimate  radial  dependence  of  relevant  ray 
parameters.  The  measured  rays  are  marked  in  the  mass-images  shown  in  the  left  column 
of  Fig.  [1]  (hereinafter  denoted  as  rays  a- e).  In  the  second  column  of  Fig.  [J] the  ray  widths, 
expressed  in  units  of  the  solar  radius  Tq,  are  presented  as  a  function  of  the  radial  distance, 
D(R),  where  R  =  t/vq.  The  ray  in  Fig.  QJi  is  characterized  by  “multi-ray”  structure,  so  we 
measured  only  the  central  feature,  which  was  the  most  prominent  element  of  the  structure. 

Inspecting  the  graphs  one  finds  that  in  the  distance  range  R  =  2-2.5  the  widths  span 
between  0.1  and  0.4  r®  (mean  0.26).  In  the  range  R  =  5-6  the  widths  increase  to  0.4- 
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Fig.  4.  The  five  measured  post-CME  rays  (indicated  by  arrows  in  the  full-resolution 
LASCO-C2  images  shown  in  the  first  column):  a)-c)  9  January  2002  at  00:06,  06:06,  and 
12:06  UT,  respectively;  d)  18  November  2003  at  12:50  UT;  e)  26  June  2005  at  06:30  UT. 
In  the  second,  third,  and  fourth  column  we  show  the  ray  width,  column-mass  excess,  and 
the  density  excess,  as  a  function  of  the  radial  distance.  The  widths  and  radial  distances 
are  expressed  in  units  of  the  solar  radius.  The  least-square  fits  are  shown  in  the  insets, 
together  with  the  correlation  coefficient  C. 
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0.9  rn  (mean  0.7).  The  widths  we  found  at  R  ~  2  are  more  than  twice  larger  than  found 


by 


Webb  et  al 


(120031)  who  analyzed  the  coronagraph  data  on  board  the  Solar  Maximum 


Mission  (SMM):  converting  their  values  of  the  width  expressed  in  degrees,  into  units  of 
the  solar  radius  at  the  corresponding  heights,  we  find  that  the  mean  width  reported  by 
Webb  et  al,  (2003j,  see  Tabs.  2  and  3  therein)  is  0.09  rQ  in  this  height  range.  The  difference 


may  be  related  to  the  lower  sensitivity  of  SMM. 


3.3.  The  ray  density 


In  the  third  column  of  Fig.  [4]  we  show  the  radial  dependencies  of  column-mass  excesses 
5m(R ),  expressed  in  g/pixel.  Inspecting  the  graphs  5m(R)  we  find  that  the  slopes  of  the 
fitted  power-law  functions  range  from  —2.1  to  —3.2.  The  average  value  adds  up  to  — 2.6±0.5. 

The  column-mass  excess  5m  was  converted  to  the  number-density  excess  5n  by  assuming 
that,  due  to  the  spherical  geometry,  the  column  length  increases  proportionally  with  the 
radial  distance,  A  oc  R.  In  particular,  we  used  Ao  =  100  Mm  at  R  =  Rq  =  2.16  (the  lowest 
height  of  measurements).  The  outcome  is  presented  in  log- log  graphs  in  the  fourth  column 
of  Fig.  |4l  together  with  the  power  law  fits.  Note  that  some  other  choice  of  Ao  would  shift 
vertically  the  5n(R)  dependence,  without  changing  the  slope  in  log-log  graph.  The  slope 
would  change  only  if  some  other  radial  dependence  A (R)  would  be  applied,  e.g.,  the  angle 
9  might  be  height-dependent. 

Inspecting  the  graphs  5n(R)  in  Fig.  [4]  we  find  that  the  slopes  of  the  fitted  power-law 
functions  range  from  —3.1  to  —4.2,  the  average  value  adding  up  —3.6  ±  0.5.  [We  pay 
special  attention  to  the  slope  of  the  fit,  since  it  does  not  depend  on  the  presumed  value 
of  Ao,  so  it  can  be  directly  compared  with  the  model  results  (Sect.  4.1).]  In  this  respect 
it  is  important  to  note  that  in  the  case  of  the  smallest  slope  in  the  sample,  the  ray  width 
becomes  approximately  constant  beyond  R  ss  3,  indicating  a  possible  decrease  of  9  with 
the  height,  i.e. ,  the  increase  of  A.  This  would  imply  that  in  this  case  the  densities  are 
overestimated  at  larger  heights,  since  5n  oc  5m/ A,  i.e.,  the  true  slope  of  5n(R)  is  probably 
somewhat  steeper. 


Since  the  evaluation  of  the  density  excess  depends  on  the  choice  of  the  column  length 
A,  in  Fig.  [5]  we  present  an  estimate  of  the  number  of  electrons  per  unit  length  of  the  ray, 
which  is  a  parameter  independent  of  A.  If  we  approximate  the  mass-difference  profiles  across 
the  ray  by  triangular  profiles,  the  mass-excess  per  unit  length  of  the  profile  amounts  to 
M  =  D5m/ 2,  from  which  we  evaluate  the  electron-number  excess  per  unit  length  of  the 
ray,  Ne. 


Inspecting  Fig. [5ji,  where  Ne(R)  is  presented  separately  for  all  five  rays  from  Fig. 01  we 
find  that  Ne(R)  either  decreases  (rays  [a],  [b],  and  [d])  or  stays  roughly  constant  (rays  [c] 
and  [e] ).  When  averaged,  the  data  from  Fig.  [5jo  show  a  dependence  similar  to  that  obtained 
bvlCiaravella  fc  Raymond!  (|2008).  who  found  that  after  a  decrease  at  low  heights,  the  value 
of  Ne  becomes  approximately  constant.  Comparing  their  Fig.  8  with  our  Fig.  [5}o  we  find 
that  the  numbers  and  the  trends  are  similar.  However,  the  difference  is  the  radial  distance 
at  which  Ne(R)  becomes  constant:  according  to  results  by  Ciaravella  fe  Raymond!  (2008) 
this  happens  at  R  «  2,  whereas  in  our  case  the  transition  occurs  at  R  sa  4.  In  Fig.  [5Jd  we 
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Fig.  5.  Number  of  electrons  per  unit  length  of  the  ray,  Ne(R):  a)  presented  separately  for 
the  five  rays  shown  in  Fig. 0] (measurements  are  labeled  in  the  legend,  following  the  sequence 
in  Fig.  0]);  b)  averaged  over  the  five  rays  (error  bars  represent  standard  deviations),  shown 
together  with  the  power-law  fit. 


present  also  the  power-law  fit,  showing  that  the  overall  trend  could  be  described  as  well  by 
Ne  (x  R-1. 


3.4.  Coronal  flows 


Coronal  regions  beneath  CMEs  are  generally  characterized  by  highly  dynamical  intricate 
structures,  revealing  complex  magnetoplasma  flows  and  waves.  Consequently,  it  is  difficult 
to  find  persistent  systematic  flow  patterns,  especially  in  periods  of  increased  solar  activity, 
when  disturbances  from  other  CMEs  affect  the  coronal  region  of  interest.  This  also  holds 
for  flow  patterns  associated  with  post-CME  rays,  which  show  permanent  morphological 
changes  and  are  strongly  affected  by  disturbances  from  distant  CMEs. 

Yet,  in  certain  situations  some  characteristic  flows  could  be  identified.  Most  frequently, 
outward-moving  inhomogeneities  along,  the  ray  are  observed,  usually  having  velocities  of 
several  hundred  kms"1  (e.g., 


Ko  et  al 


2003 


Lin  et  al 


20051).  Such  motions  are  often  in¬ 


terpreted  in  terms  of  the  reconnection  outflow,  which  may  or  may  not  be  characterized  by 
the  Alfvenic  speed  (for  a  discussion  we  refer  to  Barta  et  ah  2008,  and  references  therein). 


In  this  respect  we  note  that  sometimes  signatures  of  shrinking  loops  are  also  observed  in 


the  wake  of  CMEs  (e.g..lSheelev  &  Wang 


2002 


Sheelev  et  al 


2004 


Sheelev  fe  Wang 


2007 


and  references  therein),  consistent  with  decelerated  reconnection  downflows. 
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Fig.  6.  a)  LASCO-C2  rimning-difference  images  of  the  west-limb  ray  of  8  January  2002 
(13:54  UT),  revealing  the  inflow  into  the  ray.  The  slice  used  in  the  stack-plot  shown  in  b) 
is  outlined  by  white  rectangle  centered  at  x  =  +3.34  rQ.  b)  Stack-plot  composed  of  slices 
taken  from  successive  running-difference  images;  cc-axis  represents  the  ordinal  number  of 
the  LASCO  image  for  8  January  2002. 


Compared  to  reconnection  outflows,  detecting  signatures  of  reconnection  inflows  is  much 
more  difficult.  A  possible  example  was  reported  by  Yokovama  et,  al  (2001),  based  on  the 

(20051.  who 


Lin  et  al 


observations  by  the  EIT/SoHO.  Another  example  was  presented  by 
employed  data  from  EIT  and  UltraViolet  Coronagraph  Spectrometer  (UVCS)  on  board 
SoHO.  Recently,  Bemporad  et  al.1  (2008j)  reported  a  “side-reconnection”  in  the  aftermath  of 
a  CME,  induced  by  the  CME  expansion:  from  the  detection  of  two  converging  reconnecting 
features  at  1.7  rQ  (probably  the  CME  flank  and  the  streamer  boundary)  the  authors 


Yokqyama  et  al 


inferred  an  inflow  speed  of  3-4  km/s,  close  to  the  5  km/s  derived  by 
020011). 

In  Fig.  [G]  we  present  a  stack-plot  showing  an  example  of  the  reconnection-inflow  pattern 
observed  above  the  west  limb  on  8  January  2002  in  the  LASCO-C2  field-of-view.  Note  that 
this  post-CME  feature  is  not  the  ray  of  8-9  January  2002  presented  in  Figs.  [3]  and [4]  which 
was  on  the  opposite  side  of  the  solar  disc.  This  west-limb  event  was  not  observed  by  UVCS, 
i.e. ,  we  do  not  know  if  it  was  characterized  by  the  Fexvm  emission,  which  is  one  of  criteria 
for  identifying  post-CME  current  sheets.  It  should  be  stressed  that  we  could  not  find  a 
clear/measurable  signature  of  reconnection  inflow  in  the  three  post-CME  rays  analyzed 
in  Sects.  3.2.  and  3.3.  In  these  events  no  suitable  inhomogeneities  could  be  identified  to 
provide  tracing  of  plasma  motion  in  the  inflow  region  (to  be  discussed  in  Sect.  5). 

The  west-limb  ray  presented  in  Fig.  [Gj  formed  after  a  faint  ejection  first  observed  in 
LASCO-C2  at  02:54  UT  on  8  January  2002  (not  listed  in  the  LASCO  CME  catalog).  Each 
stack  represents  the  cut  along  a  line  parallel  to  the  y-axis  of  a  given  LASCO-C2  512  x  512- 
pixel  running-difference  image,  showing  the  pixel-intensities  in  the  range  x  =  3.30  3.38  rQ. 
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Note  that  the  ray  was  oriented  very  close  to  the  x-axis  of  LASCO  images,  i.e. ,  the  stacks 
are  roughly  perpendicular  to  the  ray  direction.  The  stack-plot  contains  all  images  recorded 
on  8  January  2002  (denoted  by  ordinal  number  at  the  x— axis  of  the  stack-plot;  the  time 
cadence  is  on  average  24  min). 

The  black/white  stripes  converging  towards  the  ray  location,  located  at  y  ~  400"  ~ 
0.41rQ,  clearly  show  inflows  into  the  ray.  The  corresponding  velocities  in  the  left  part  of 
the  stack-plot  are  found  in  the  range  15-25  kms-1  (mean  19  ±  4  kms-1).  Such  velocities 
are  several  times  higher  than  5  kms-1  reported  by  Yokoyama  et  al.  ((20011) .  After  the  slice 
No.  50  (18:54  UT),  the  stack-plot  reveals  effects  of  a  large-scale  perturbation  from  another 
CME,  which  occurred  above  the  E-limb.  Note  that  the  “push”  caused  by  the  perturbation 
increased  the  inflow  speed,  now  ranging  from  25  to  30  kms-1  (mean  27  ±  5  kms-1). 

Finally,  it  should  be  noted  that  we  also  observed  intermittent  downflows  in  the  form 
of  shrinking  loops  (similar  to  that  reported  by  ISheelev  et  al.  120041).  starting  from  R  ~ 


3.  At  larger  heights,  outflows  could  be  recognized,  sometimes  showing  a  pattern  like  in 
disconnection  events  (e.g., 


Webb  fc  Cliver 


1995: 


Simnett  et  al 


19971 


Wang  et  al.  1999). 


This  implies  that  the  diffusion  region  in  this  event  was  located  around  R 
higher  than  in  the  events  shown  in  Fig.  0 


3,  i.e.,  much 


4.  Comparison  with  the  model  and  previous  studies 

4.1.  Comparison  with  the  model  results 

In  Fig.  0i  we  present  a  comparison  of  measurements  presented  in  4th  column  of  Fig.  [4] 
with  the  model-dependencies  Sn(R)  derived  in  the  Appendix.  In  particular,  we  present  the 
model  results  based  on  the  isothermal  (T  =  106  K)  [Parker  (|1958l)  solar  wind  model  (see  also 


Mann  et  al 


1999al).  For  the  magnetic  field  B(R)  we  take  the  empirical  coronal  magnetic 


field  scaling  established  by  Bulk  &  McLean  (1978).  Note  that  data  [a],  [b],  [d],  and  [e]  are 
slightly  shifted  in  JCcoordinate  (symmetrically  with  respect  to  [c])  to  avoid  overlapping. 
Figure  0i  clearly  shows  that  the  observations  fit  well  to  the  model  curves  (as  stated  in 
Sect.  3.3  the  choice  of  the  column  length  Ao  only  shifts  the  values  up  or  down,  but  does 
not  change  the  slope).  Comparing  the  trend  of  the  data  we  see  that  the  slope  corresponds 
much  better  to  the  model  slopes  for  the  CSs  than  for  the  ambient  corona.  Furthermore  one 
finds  that  the  CS  densities  are  more  than  one  order  of  magnitude  larger  than  that  in  the 
ambient  corona. 

In  Fig.  0>  we  present  the  measured  average  density  excess  5n(R)  (circles  with  error  bars 
represent  the  mean  values  from  Fig.  0i  with  the  associated  standard  deviations)  and  the 
corresponding  power-law  fit,  compared  with  the  CS  model  dependencies.  The  observational 
data  show  the  dependence  i?-3'6,  whereas  in  the  same  height  range  the  CS  model  data 
show  ~  R~30  -  R~3  2  and  the  quiet-corona  model  behaves  as  R~ 4  3  (see  Appendix). 


4.2.  Ray  densities  from  previous  studies 

Post-CME  rays  were  studied  in  a  number  of  papers,  employing  various  techniques  to  esti¬ 
mate  their  density.  In  Table0we  present  an  overview  of  the  spectrographic  results  obtained 
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Fig.  7.  a)  Comparison  of  the  density-excess  measurements  presented  in  the  right  hand 
column  of  Fig.@]with  the  model  results  (thin  black  lines  labeled  by  the  heliocentric  distance 
of  the  diffusion  region,  i?x).  Measurements  are  labeled  in  the  legend,  following  the  sequence 


in  Fig.  [U  The  thick-gray  line  shows  the  density  model  of  a  quiet  106  K  corona  ([Mann  et  al 


1999a).  b)  The  mean  values  of  the  data  presented  in  (a),  shown  by  circles  with  error  bars  and 


the  corresponding  power-law  fit  (written  in  the  inset),  compared  with  the  model  results  and 
the  values  from  previous  studies  reviewed  in  Sect.  4.2  (diamond  -  . 


Ciaravella  et  al 


(120021);  squares  - 


Ko  et  al 


Ko  et  al 


2003);  triangle 


(20031);  cross  -  Bemporad  et  al.  (20061)). 


from  UltraYiolet  (UV)  spectra  of  five  post-CME  rays  (dates  are  given  in  the  header  row 


reported  by 


Ciaravella  et  al 


(20021) 


Ko  et  al. 


(2003b 


Raymond  et  alJ  ( 20031) . 


Lee  et  al 


(2006).  Bemporad  et  ah!  ((2006),  andlciaravella  fe  Ravmondl  (koosl).  In  the  last  column  we 
also  present  the  unpublished  results  by  Schettino  et  al.  (2008,  to  be  submitted). 

All  rays  listed  in  Table  |T]  were  characterized  by  the  presence  of  the  narrow  Fe  XVIII 
emission.  UV  spectra  of  rays  have  been  obtained  with  the  UltraViolet  Coronagraph 
Spectrometer  (UVCS;  Kohl  et  al,  19951)  on  board  SoHO.  The  spectrometer  observes 


through  a  narrow  entrance  slit  42'  long  and  up  to  84"  wide.  The  slit  can  be  located  at 
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Table  1.  Spectrographic  results  on  the  UV  current  sheets. 


1998  Mar  231 

2002  Jan  82 

2002  Apr  213 

2002  Nov  264 

2003  Nov  45 

2003  June  26 

CME  lst-LASCO 

09:33:36 

17:54:05 

01:27:20 

17:06:21 

19:54:05 

08:54:05 

CME  speed  (kms-1) 

403 

1794 

2393 

479 

2657 

980 

Flare 

N 

Ybl* 

X  1.5 

N 

X  30.8 

M3.9 

AR  No. 

— 

9782/85 

9906 

— 

10486 

10365 

LASCO  CS 

N 

2  days 

N 

N 

~  20  h 

6.5  h 

CS  obs.  start 

16:00:30 

17:48/Jan  10 

00:45:34 

18:39:15 

20:03:50 

09:40:44 

CS  obs.  height 

1.51  (-1.38) 

1.53 

1.62 

1.61 

1.69 

1.68 

CS  FWHM  width  (rQ) 

0.1 

0.2 

0.3 

0.4 

0. 2-0.1 

0.3 

CS  full  width  (r©) 

0.4 

0.4 

0.6 

0.8 

0.5-0. 2 

0.35 

CS  LoS  depth  (r©) 

0.06 

0.2 

— 

0.5 

0.1 

0.1 

CS  obs.  duration 

20  h 

10  h 

14  min 

2.3  days 

~  17  h 

6.3  h 

CS  temp.  (106  K) 

5 

3-4 

5 

8-3 

8-4 

7 

CS  density(107cm-3) 

5-10 

<  4 

— 

6.5 -7.5 

7-10 

1 

1  - 
2  - 

3  - 

4  - 

5  - 


Ko  et  al.  (200. 


Ciaravella  et  al. 


3); 


12002); 


Ra^mond et-^d 


(2003) 


Lee  et  al. 


(20061); 


Bemporad  et  al.  (2006); 


Ciaravella  fc  Raymond!  (|2008l) 


6  -  Schettino  et  al.  (2008,  to  be  submitted) 
*  Ybl  =  Yes  (behind  the  limb) 


any  polar  angle  and  at  heliocentric  distances  from  1.5  up  to  10  rg.  UVCS  can  image  spec¬ 
tra  of  the  solar  corona  in  the  range  945-1270  A  (473-635  A  in  the  second  order). 

The  first  five  rows  of  Table  1  specify  the  associated  CME/flare  events:  the  first  appear¬ 
ance  of  the  CME  in  the  LASCO-C2,  the  CME  mean  speed,  the  soft  X-ray  importance  of 
the  flare  if  observed,  the  active  region  label,  and  duration  of  the  LASCO-ray  if  observed, 
respectively. 

The  next  eight  rows  concern  the  ray  characteristics  as  detected  in  the  UV  spectra:  the 
time  of  the  first  detection  of  the  ray  by  UVCS,  the  height  at  which  the  spectrograph  slit 
intersects  the  ray,  the  ray  width  along  the  UVCS  slit  expressed  in  units  of  Tq  (both  the 
full- width  at  half-maximum  and  the  full  width  are  presented),  the  estimate  of  the  LoS 
depth,  duration  of  the  ray  UVCS  observations,  and  finally,  in  last  two  rows  we  present  the 
ray  temperature  and  density. 

The  ray  UVCS  densities  are  presented  in  Fig.  [7b,  together  with  the  density-excess 
measurements  presented  in  Sects.  3.3  and  4.1.  The  estimate  based  on  the  UVCS  spectra 


reported  by 


Ko  et  al 


(2003)  is  indicated  by  the  diamond  symbol.  The  column- length  was 


assumed  therein  to  be  140  Mm  which  was  taken  to  be  the  same  as  the  FWHM  of  the  spatial 
profile  of  the  FeXVIII  A974  emission  across  the  post-CME  rays.  Such  a  column-length  is 
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about  2  times  larger  than  the  value  we  would  get  for  this  height  by  the  scaling  that  was 
applied  in  Sect.  3.  Thus,  to  adjust  the  values  to  our  Ao  the  density  has  to  be  multiplied  by 
a  factor  of  2  (indicated  by  vertical  bar  attached  to  the  diamond  symbol  in  Fig.  [TJo)  - 


The  density  reported  by 


Ciaravella  et  al 


is  shown  by  the  black  triangle. 


(20021.  also  estimated  from  the  UVCS  data, 


Ciaravella  et  al 


(2002)  estimated  the  CS  density  at  R  =  1.5 


to  5-  10  x  107  cm-3  assuming  the  column  length  of  40  Mm.  After  adjusting  to  the  column- 
length  equivalent  to  that  we  used  in  Sect.  3  (required  increase  is  around  20  %),  we  find  the 
equivalent  density  around  6  x  10'  cm~3. 

The  density  inferred  by  Bemporad  et  al.  (20011)  from  the  spectral  observations  of  the 
face-on  current  sheet  of  26  November  2002  is  shown  by  the  black  cross.  The  attached 
vertical  bar  indicates  that  the  density  could  be  lower  if  a  larger  current  sheet  thickness 
would  be  assumed.  The  gray  plus  symbol  represents  the  estimated  density  of  the  ambient 
coronal  plasma. 

In  Fig.  [7):)  we  also  show  (black  squares)  the  density  excess  data  based  on  the  ray  and 


coronal  densities  estimated  by 


Ko  et  al. 


(2003)  employing  the  LASCO-C2  observations  and 


using  the  Thompson  scattering  function  ([Billings  1966j) .  In  estimating  the  ray  density  at 
(J2003I)  assumed  the  column-length  of  ~  250  Mm,  which  is  about  1.2 


R  =  4.4 


Ko  et  al 


times  larger  value  than  we  used  at  the  same  height.  Thus,  adjusting  their  results  to  our 
column-length  would  result  in  1.2  times  lower  density. 

Finally,  in  Fig.  [T]>  the  results  based  on  the  polarization-brightness  measurements  per¬ 


formed  by 


Poletto  et  al 


(12008)  are  presented.  The  thick  black  line  represents  the  density 


excess  Sn  re-evaluated  from  Fig.  9  of 


Poletto  et  al 


(20081)  by  applying  column-lengths  con¬ 


sistent  with  those  applied  in  Sect.  3.  The  ray  data  drawn  by  the  thick  line  fit  very  closely 
to  our  measurements,  whereas  the  “quiet  corona”  data  (thin  gray  dashed  line)  are  very 
close  to  the  model  densities. 

The  power-law  fit  to  all  data  (adjusted  to  the  same  Ao)  shown  in  Fig.  [T}o  has  the  slope 
R~3-3  (without  the  Ao  adjustment  would  add  up  to  R~3  0).  This  is  quite  close  to  the 
model  slopes  found  in  the  radial  distance  range  1.5-7  r^,  characterized  by  the  power-law 
exponents  3.0 -3.2  for  l?x  =  1.1  1.5  and  the 


Dulk  fe  McLeanl  ( 1978)  magnetic  field. 


5.  Discussion  and  conclusion 

The  empirical  characteristics  of  post-CME  rays  we  summarize  as  follows: 

1.  rays  often  show  activity  in  the  form  of  outflowing  features,  morphological  changes,  and 
changes  of  the  inclination;  sometimes  reconnection  inflows  are  observed; 

2.  the  width  of  rays  increases  with  height,  from  ~0.1  0.3  r0  at  R  ~  2  to  ~ 0.4-0. 8  r0 
at  R  ~  6. 

3.  densities  found  in  rays  are  at  least  several  times  (up  to  more  than  one  order  of  magni¬ 
tude)  larger  than  in  the  ambient  corona  in  the  considered  height  range; 

4.  coronal  regions  surrounding  rays  are  depleted; 

5.  on  average,  the  number  of  electrons  per  unit  length  of  the  ray,  Ne,  first  decreases  with 
the  height,  and  then,  at  larger  heights  becomes  approximately  constant; 
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6.  temperature  of  rays  in  the  range  R  ~  1.5 -1.7  spans  from  3  to  8  MK,  i.e. ,  it  is  several 
times  larger  than  in  the  “normal”  corona,  with  a  tendency  to  decrease  in  time. 

The  nature  of  the  coronagraphic  white-light  ray  activity  is  twofold.  The  formation  of 
fine  structure  elements,  mostly  blob-like  features,  and  their  outward  motion  at  speeds  in 
the  range  100-  1000  kms-1,  seems  to  be  the  most  common  internal  activity.  It  could  be 


Bart  a  et  al. 

2008; 

■  07 - -  ~ 

Rilev  et  al. 

2007 

ejection  of  plasmoids  (e.g., 

On  the  other  hand,  changes  of  the  ray  inclination,  as  well  as  morphological  changes  that 
could  be  attributed  to  changes  of  the  ray  geometry,  seem  to  be  caused  by  the  large-scale 
magnetic  field  evolution  in  the  wake  of  the  CME.  In  some  cases  this  type  of  activity  is 
caused  by  perturbations  coming  from  distant  eruptions,  most  likely  the  large-amplitude 
waves  or  shocks. 

As  stated  in  item  2,  the  width  increases  on  average  by  about  0.4  rQ  over  the  distance  of 
4 Tq,  i.e.,  the  width  is  roughly  proportional  to  the  radial  distance.  Following  our  interpre¬ 
tation  of  rays  in  terms  of  the  reconnecting  current  sheet,  the  ray  boundaries  should  outline 
the  slow  mode  shocks.  Given  the  model  presented  in  the  Appendix,  in  particular  Fig.  [HI  it 
can  be  concluded  that  the  slow  mode  shocks  are  oriented  very  close  to  the  radial  direction, 
i.e.,  that  the  angle  denoted  in  Fig.  [8]  as  </>,  is  very  small.  Since  the  angle  (j)  (expressed  in 
radians)  equals  to  the  inflow  Alfven  Mach  number  reduced  by  the  factor  712/711  that  rep¬ 
resents  the  density  jump  at  the  slow  mode  shock  (the  value  of  712/711  depends  on  f3,  and 


should  be  around  2;  see  Vrsnak  &  Skgnder  2005),  one  finds  that  the  inflow  Mach  number 
is  low,  most  likely  in  the  order  of  0.01  or  less.  [For  larger  values,  (j>  would  be  large  enough 
to  make  the  “super-radial”  widening  of  the  ray  measurable.]  This  may  be  the  reason  why 
inflows  like  those  shown  in  Fig.  [5]  are  rarely  observed. 

Note  that  the  situation  shown  in  Fig.  [G]  is  different  from  the  previously  described  “sta¬ 
tionary”  pattern.  In  this  case  we  see  distinct  elongated  features  inflowing  from  both  sides 
towards  the  ray  axis  of  symmetry.  Thus,  these  features  cannot  be  interpreted  as  CS  bound¬ 
aries,  i.e.,  signatures  of  the  slow  mode  shocks,  but  rather  they  represent  coronal  density 
inhomogeneities  aligned  with  the  magnetic  fieldlines.  These  features  are  inclined  from  the 
radial  direction,  i.e.,  show  the  “super-radial”  orientation.  When  observed  in  pairs,  they 
form  a  V-pattern  roughly  symmetric  with  respect  to  the  ray  axis.  Thus,  these  substruc¬ 
tures  most  likely  outline  the  magnetic  field  inflowing  into  the  current  sheet.  In  such  a 
situation  the  angle  (expressed  in  radians)  between  the  fieldline  and  the  axis  of  symmetry  is 
equal  to  the  inflow  Alfven  Mach  number  Ma  (e.g.,  Vrsnak  &  Skender  2005).  In  the  event 


shown  in  Fig.  [6l  we  estimate  this  angle  to  a  few  degrees,  so  we  can  take  Ma  ~0.02  0.05. 
Bearing  in  mind  that  the  inflow  speed  was  estimated  to  ~  20  kms-1,  this  would  correspond 
to  the  Alfven  speed  of  400  1000  kms-1,  which  seems  reasonable  for  this  height  range  (see, 


e-g-, 


Vrsnak  et  al 


200411. 


The  main  physical  characteristics  of  post-CME  rays  are  their  mass/density  excess  and 
increased  temperature  (Table  I).  This,  as  well  as  the  morphology  and  flows,  can  be  explained 
in  terms  of  the  reconnection  outflow-jet,  being  the  structural  element  of  the  vertical  current 
sheet  that  forms  in  the  wake  of  a  CME.  The  measured  values  of  the  CS  density  and  the 
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temperature  provide  also  an  estimate  of  the  ambient  coronal  magnetic  field,  since  the 
external  magnetic  pressure  is  roughly  equal  to  the  gas  pressure  in  the  outflow  region  (see, 
e.g.,  Appendix  inlAurass  et  al.ll2002ll.  Using  the  values  ncs  =  5-10  xlO7  cm-3  and  Tcs  = 
5-8  MK,  we  find  for  the  ambient  magnetic  field  Bc  ~  0.9-  1.7  gauss.  Such  a  magnetic  field 


Dulk  fc  McLean 


is  consistent  with  the  empirical  scaling  B  =  0.5 (R—  l)-1'5  established  by 
(1978|),  which  gives  B  ~  0.85-1.41  gauss  for  the  radial  distance  range  R  =  1.5 -1.7. 
Similarly,  using  the  relationship  between  the  temperature  jump  at  SMSs  and  the  external 


Aurass  et  al 


20021. 


plasma-to-magnetic  pressure  ratio,  T2/Ti  =  1  +  0.4//?  (see  Eq.  12  in 
one  can  estimate  the  value  of  /?  in  the  inflow  region.  Applying  T2/T\  <  Tcs/Tc  =  2-8,  we 
find  /?  >  0.4-  0.06.  Such  values  are  roughly  consistent  with  those  expected  for  the  active 
region  corona  in  this  height  range  (IGarvil200llh 

The  results  of  the  model,  based  on  the  quasi-stationary  Petschek  (|1964f)  reconnection 
regime,  are  consistent  with  the  observed  density  excesses.  At  the  standing  slow  mode  shocks, 
formed  by  the  CS  inflow,  the  plasma  is  heated,  compressed,  and  accelerated/deflected  to 
form  the  upward  directed  reconnection  jet.  In  this  way  the  dense  plasma  from  low  corona 
is  transported  to  larger  heights,  causing  the  observed  density  excess.  Model  results  show 
that  the  density  can  be  increased  by  more  than  one  order  of  magnitude,  which  is  consistent 
with  observations.  However,  it  should  be  noted  that  the  observed  values  depend  on  the 
assumed  value  of  the  LoS  CS-depth  Ao-  Our  results  indicate  that  Ao  is  in  the  order  of  100 
Mm,  and  that  the  plane  of  the  observed  rays  is  inclined  at  small  angle  with  respect  to  the 
LoS. 

The  parameter  that  does  not  depend  on  the  LoS  CS-depth,  or  the  CS  orientation,  is 
the  number  of  electrons  per  unit  length  of  the  CS,  which  we  denoted  in  Sect.  3.3  as  Ne. 
The  measurements  show  that  Ne(R)  either  decreases,  or  is  approximately  constant.  On 
average,  Ne(R)  first  decreases,  to  become  approximately  constant  at  Rc  ~  4.  Given  our 
observations,  as  well  as  the  observations  by  [Ciaravella  fe  Raymond!  (l2008ll).  it  seems  that 
Rc  is  different  in  different  events. 

Bearing  in  mind  the  equation  of  continuity,  there  are  two  effects  that  determine  the 
radial  dependence  of  Ne.  One  is  the  radial  dependence  of  the  reconnection-outflow  velocity, 
and  another  one  is  the  contribution  of  the  inflowing  plasma.  For  example,  if  the  outflow 
speed  would  be  constant,  the  value  of  Ne  would  be  increasing  monotonously  with  the  height 
due  to  the  cumulative  supply  of  the  plasma  through  the  SMSs  (more  and  more  inflowing 
plasma  joining  to  the  outflow).  On  the  other  hand,  if  the  inflow  contribution  is  negligible 
(as  it  is  at  large  heights),  then  a  decelerated  outflow  would  be  associated  with  increasing 
Ne(R ),  whereas  the  accelerated  outflow  would  lead  to  a  decreasing  Ne(R). 

Inspecting  outcomes  for  various  model  inputs,  we  found  that  regardless  on  the  model 
details,  in  the  vicinity  of  the  diffusion  region  Ne(R)  increases  steeply  due  to  a  strong  effect 
of  the  plasma  inflow.  At  larger  heights,  beyond  ~O.5r0  above  the  diffusion  region,  the 
behavior  of  Ne(R)  becomes  dependent  on  the  model  input.  A  match  with  the  observations 
(decreasing  Ne(R)  followed  by  Ne  ~  const.)  is  found  only  for  cases  where  the  Alfven 
velocity  and  the  solar  wind  speed  increase  in  the  given  height  range,  and  only  if  the  diffusion 
region  is  below  R  ~  1.4.  Given  that  the  power-law  fit  presented  in  Fig.  [5})  quite  well 
describes  the  Ne(R)  dependence  as  Ne  oc  f?_1,  would  imply  that  the  outflow  velocity 
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increases  approximately  as  vout  oc  R.  However,  note  that  the  last  three  data-points  in 
Fig.  [5)3  might  indicate  that  Ne  starts  to  increase  gradually  beyond  R  ~  5.  Such  a  behavior 
is  found  in  all  of  the  considered  model  options,  i.e. ,  after  the  solar  wind  speed  becomes 
approximately  constant  and  the  Alfven  velocity  starts  decreasing,  the  value  of  Ne  gradually 
increases  with  the  radial  distance. 

In  this  respect,  it  should  be  stressed  that  the  model  results  generally  depend  on  the 
applied  coronal  density  and  magnetic  field  model.  On  the  other  hand,  the  Ne(R)  measure¬ 
ments  are  burdened  by  large  errors  in  estimating  the  ray  width  D.  So,  since  the  comparison 
of  the  observed  and  the  calculated  Ne(R)  might  lead  to  ambiguous  conclusions,  it  turns 
out  that  the  most  reliable  model-parameter  to  be  compared  with  observations  is  the  slope 
of  the  ncs(R)  dependence,  which  should  be  less  steep  than  the  slope  of  the  ambient  nc(R), 
regardless  on  the  model  details.  Comparing  the  calculated  and  observed  slopes  we  find  a 
very  good  correspondence. 

Finally,  the  model  predicts  a  decrease  of  the  density  excess  in  time,  due  to  the  rise  of 
the  diffusion  region  height.  Bearing  in  mind  that  this  rise  is  slow,  it  is  to  be  expected  that 
the  decrease  of  the  density  excess  should  be  slow  too.  Indeed,  we  did  not  detect  a  decrease 
of  the  ray  density  in  12  hour  interval  analyzed  in  the  ray  of  9  January  2002.  However,  note 
that  a  decrease  of  the  density  was  reported  by  ICiaravella  &  Raymond!  (]2008[)  in  the  event 
of  4  November  2003.  From  their  Table  1  we  find  that  the  density  decreased  by  ~  40-50% 
in  about  2  hours.  According  to  our  model,  such  a  decrease  would  correspond  to  an  increase 
of  the  diffusion  region  height  for  A R  ~  0.1,  which  would  correspond  to  the  rise  speed  in 
the  order  of  10  kms-1.  On  the  other  hand,  given  the  accuracy  of  our  density  estimates 
(~  50%),  we  infer  that  in  the  event  of  9  January  2002  the  rise  speed  of  the  diffusion  region 
was  less  than  2  kms-1. 


Appendix 

The  overall  geometry  of  the  post-CME  ray,  assumed  to  be  a  signature  of  the  bifurcated 
reconnecting  current  sheet,  is  shown  in  Fig.  [TJ  The  diffusion  region,  where  the  magnetic 
field  lines  reconnect,  is  located  at  the  radial  distance  Rx.  The  oppositely  directed  field 
lines  merge  at  the  velocity  Vin(R),  bringing  into  the  reconnection  outflow  jet  (the  ray)  the 
coronal  plasma  of  the  density  nc(R).  We  assume  that  the  merging  velocity  is  faster  than  the 
local  slow- mode  speed,  so  the  slow  mode  shocks  (SMSs)  form  in  between  the  two  inflows, 
bounding  the  reconnection  outflow  like  in  the  Petschek  1964|)  reconnection  model. 

Given  the  spherical  geometry  of  the  corona,  the  flow  velocity  vector  in  the  inflow  region 
can  be  represented  as  a  superposition  of  the  radial  component  (corresponding  to  the  solar 
wind  speed  vsw)  and  the  “horizontal”  component  t>in)  (direction  of  the  unit  vector  $  shown 
in  the  lower-right  corner  of  Fig.  [8]),  as  depicted  in  Fig.  [[):>■ 

In  the  stationary  state  the  relation  V  x  E  =  0  must  be  satisfied,  to  provide  <9B /dt  =  0.  In 
a  spherical  coordinate  system  (see  the  lower-right  corner  of  Fig.  this  implies  d{rE,f,)  /  dr  = 
0,  i.e.,  Rv\„B  =  const.  For  the  ambient  coronal  magnetic  field  B(R)  we  take  the  empirical 


]see  also 


Vrsnak  et  al 


model  by  Dulk=&^cLean  197811  and  the  simple  analytical  model  by 


Mann  et  al 


11999b) 


(2002)],  which  then  defines  i’in(R). 
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We  denote  the  angle  between  the  radial  direction  and  the  slow  mode  shock  (Fig.  [HI)  as 
4>{R).  Taking  into  account  the  continuity  equation  across  the  SMS 


vinpcAR  =  tW/OSMsAi?.0  , 


(1) 


where  r^ut  ~  va  +  vsw  ( Skender  et  al 


2003 


Vrsnak  fc  Skender 


200511.  we  find 


0  =  — - j*—  >  (2) 

PSMS  VA  +  Vsw 

which  defines  the  geometry  of  the  SMSs. 

For  the  coronal  density  pc(R)  we  use  the  isothermal  solar  wind  model  by  Parker  (see 


Mann  et  al 


1999af)  using  for  the  temperature  Tc  =  1,  1.5,  and  2  MK,  which  provides  also 


the  corresponding  solar  wind  speed  vsw(R).  This,  together  with  the  previously  defined 
B(R ),  also  defines  the  ambient  coronal  Alfven  velocity  va(R )  and  plasma-to-magnetic 
pressure  ratio  (3(R),  both  governing  the  jump  relations  at  the  SMS.  Assuming  that  the 
guiding  field  is  negligible,  the  density  and  temperature  jump  at  SMS  can  be  expressed  as: 

PSMS  5(1  +  /?) 


(Aurass  et  al 


2002; 


Skender  et  al 


Pc 

Isms 

—L, 


2 

=  1 


2003; 


-5/3  ’ 
2 

5/?l 


(3) 


(4) 


Vrsnak  fc  Skender 


2005).  which  determines  the 


density  and  the  temperature  of  the  plasma  inflowing  into  the  reconnection  outflow  jet  at  a 
given  R. 

At  a  given  radial  distance  R  the  density  of  the  reconnection  outflow  is  determined  by 
the  flow  carried  from  lower  heights  and  by  the  flow  carried  across  the  SMS.  Bearing  in 
mind  the  geometry  depicted  in  Kig.  [Sj  the  continuity  equation  (mass  conservation)  can  be 
written  in  the  form: 


p°?lV%\Li+1di+1  =  pT^Udi  +  p?+1v™+1Li+1AR ,  (5) 

where  Li+ 1  =  (Li  +  Li+ 1)/2,  p1°i_1  =  pln(Ri+ 1),  and  v ^  =  vln(Ri+±)  are  the  current 
sheet  length,  the  ambient  density,  and  the  ambient  velocity  at  R  =  Ri+i  =  Ri  +  AR/2 , 
respectively.  The  term  on  the  left-hand  side  represents  the  mass  flow  through  the  upper 
surface  Ti+i  x  di+ 1,  where  L  represents  the  length  of  the  current  sheet  perpendicular  to 
the  plane  of  Fig.  [8]  The  first  term  on  the  right-hand  side  represents  the  mass  flow  through 
the  bottom  surface  Li  x  di  and  the  second  one  the  inflow  across  the  SMS.  The  flow  speeds 
vout  can  be  estimated  as  v°ut  =  va (Ri)  +  vsw(Ri)  and  =  VA(Ri+i)  +  vsw(Rl+i). 
Bearing  in  mind  the  spherical  geometry,  the  length  L  reads 

L(R)  =  L(RX)^~  ,  (6) 

Jlx 

where  Rx  is  the  radial  distance  of  the  diffusion  region  and  L(RX)  is  the  current  sheet 
extension  in  the  ^-direction  at  Rx,  i.e.,  it  represents  the  length  of  the  current  sheet  X-line. 
The  width  of  the  current  sheet  at  a  given  radial  distance  can  be  evaluated  using 

di+ 1  =  di  1~  +  ARcpi+i ,  (7) 

41  i 

with  dx  =  d(Rx)  ~  0,  i.e.,  d\  =  (j)\AR.  The  first  term  on  the  right-hand  side  of  Eq.  (JTJ)  is 
due  to  the  spherical  geometry  (hereinafter  we  abbreviate  diRi+i/Ri  =  d*+11  i.e.,  di+i  = 
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d*+1  +  A.R(f>i+i),  and  the  second  term  is  due  to  the  inclination  of  the  SMS  with  respect  to 
the  local  radial  direction  (Fig.  [5}. 

Equations  (Hhzd  provide  evaluation  of  the  density  distribution  along  the  post  CME- 
ray,  p0ut(R),  for  various  values  of  diffusion  region  height  /?y.  In  Fig  [5]  we  show  the  results 


Dulk  &  McLean  (  and 


Mann  et  al 


obtained  using  the  coronal  magnetic  field  model  by 
(11999b).  whereas  for  the  coronal  density  we  use  the  isothermal  (106  K)  solar  wind  model 


by 


Mann  et  al 


dl999aj). 


Let  us  now  consider  the  plasma  temperature  along  the  current  sheet.  In  the  first  layer 
above  the  X-line  the  temperature  equals  the  temperature  Tsms  evaluated  by  Eq.  ©  at 
the  radial  distance  R  =  R±  =  Rx  +  AR/2.  The  plasma  of  this  temperature  moves  outward 
and  enters  into  the  next  layer  through  the  area  d±  x  L i,  undergoing  adiabatic  cooling  due 
to  volume  expansion.  In  the  second  layer  it  meets  with  the  plasma  which  inflows  into  the 
current  sheet  through  the  slow  mode  shock,  where  it  was  heated  to  Tsms(-R2)-  The  same 
happens  successively  in  each  further  level. 

Since  thermal  conductivity  of  coronal  plasma  is  very  high,  the  temperature  can  be 
consider  as  uniform  across  each  layer,  i.e. ,  along  the  magnetic  field  lines  connecting  SMSs 
(magnetic  field  component  B$).  Thus,  we  can  write  for  the  temperature  of  a  given  layer: 

rridTd  +  tosmsTsms 


Tcs  = 


(8) 


md  +  tosms 

Here,  md  and  Td  are  the  mass  and  temperature  of  the  plasma  contained  in  the  volume 
defined  by  the  lines  denoted  in  Fig.  [8]  as  di  and  d*+1.  Note  that  we  have  neglected  the 
thermal  conductivity  along  the  CS  axis  of  symmetry,  since  the  heat  flow  in  this  direction  is 
reduced  due  to  the  magnetic  field  component  B$,  and  the  temperature  gradients  are  small 
due  to  large  vertical  length-scale  of  CS.  The  mass  md  is  determined  by: 


A+i 


Pi+l  +  Li+1  +  Li 


-S+l 


di 


(9) 


a  2  2  2 
whereas  tosms  and  Tsms  are  the  mass  and  temperature  of  the  plasma  contained  in  the 
“triangular”  volume  between  the  radial  and  the  slow  mode  shock.  The  temperature  Tsms 
can  be  evaluated  by  Eq.  (J4])  at  the  radial  distance  R  =  Ri+i  =  Ri  +  AR/2,  whereas  the 
mass  tosms  is  determined  by: 


mSMS  —  Psms^^ 


i+ 1 


Li+ 1  +  Li  AR(j)i+i 


(10) 


2  2 

where  pg^g  is  evaluated  by  Eq.  Q  at  the  radial  distance  R  =  Ri+i  =  Ri  +  AR/2. 

The  estimate  of  the  temperature  Tcs  in  a  given  layer  still  requires  evaluation  of  the 
temperature  Td,  where  we  have  to  take  into  account  the  adiabatic  expansion  from  the  layer 
“i”  to  “i  +  1”.  The  plasma  contained  in  the  volume  defined  by  lines  di- 1  and  df. 

Li  +  Li- 1  di  +  di- 1 


Vi  =  AR 

expands  into  the  volume: 

Li+ 1  +  Li  d*+1  +  di 


Vt+i  =  A  R 


The  time  At  can  be  expressed  as 


At  = 


+  (Li+1d*+1v°»\  -  LidiV^ 


A  R 


)A  t. 


Ji+ 1 


(ii) 


(12) 


(13) 
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Fig.  8.  Geometry  of  the  post-CME  ray  model.  Dash-dotted  lines  depict  the  axis  of  sym¬ 
metry,  shaded  area  represents  the  slow  mode  shock  (SMS)  and  dotted  lines  the  local  radial 
direction.  Local  coordinate  system  is  indicated  in  the  lower-right  corner.  For  details  see  the 
text  in  the  Appendix. 

where  =  (i>°b*  +  v°ut)/2.  The  difference  of  the  first  term  on  the  right-hand-side  of 
Eq.  (1121)  and  Vi  represents  the  expansion  in  ip  and  •&  direction,  whereas  the  second  term 
on  the  right-hand-side  of  Eq.  m  represents  the  contribution  of  the  expansion  in  the 
f-direction.  Bearing  in  mind  the  characteristics  of  the  adiabatic  expansion,  finally  we  get: 


(14) 


where  7  is  the  ratio  of  specific  heats,  for  hydrogen  plasma  being  equal  to  7  =  5/3. 

The  system  Eqs.  (PHI  can  be  solved  starting  from  the  first  layer  above  the  diffu¬ 
sion  region,  and  successively  calculating  the  values  in  every  new  layer  by  using  the  values 
obtained  for  previous  layer.  In  this  way  we  get  the  dependencies  p{R)  and  T(R)  for  the 
current  sheet  plasma. 
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